Paramagnetic Breakdown of Superconductivity in Ultrasmall Metallic Grains 
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We study the magnetic-field-induced breakdown of superconductivity in nm-scale metal grains 
having a mean electron level spacing d ~ A (bulk gap). Using a generalized variational BCS approach 
that yields good qualitative agreement with measured spectra, we argue that Pauli paramagnetism 
dominates orbital diamagnetism, as in the case of thin films in a parallel magnetic field. However, the 
first-order transition observed for the latter can be made continuous by finite size effects. The mean- 
field procedure of describing the system by a single pairing parameter A breaks down for d ~ A. 



PACS numbers: 74.20.Fg, 74.25.Ha, 74.80.Fp 

When a system of (correlated) electrons is sufficiently 
small, the electronic spectrum becomes discrete. This al- 
lows one to study the nature of electron correlations in 
unprecedented detail by analyzing the details of the spec- 
trum. It has recently become possible to measure such 
discrete spectra directly by studying electron transport 
through nm-scale metallic grains (radius r « 5nm), for 
which the mean spacing d = l/N{eF) is ~ 0.1 — 1 meV 
For Al grains the effects on the spectrum of spin- 
orbit interactions R], non-equilibrium excitations |^ and 
superconductivity y|^J|] have been investigated. 

Studying the latter is particularly interesting in grains 
with d ~ A (bulk gap) , near the lower size limit |^ of ob- 
servable superconductivity. The number of free-electron 
states that pair-correlate (those within A of Ef) is then of 
order one. Thus, even in grains in which a gap can still be 
observed j^] , pairing correlations are expected to become 
so weak that they might be destroyed by the presence 
of a single unpaired electron Q]. A direct way to probe 
this is to turn on a magnetic field, whose Zeeman energy 
favors paramagnetic states with non-zero total spin. 

In this Letter, we develop a theory for the paramag- 
netic breakdown of superconductivity in nm-scale grains. 
We exploit analogies to thin films in a parallel magnetic 
field , but explicitly take account of the discreteness of 
the spectrum. To calculate the eigenenergies En of the 
grain's lowest-lying eigenstates \n), we adopt a general- 
ized variational BCS approach that goes beyond standard 
mean-field theory by using a different pairing parameter 
A„ for each \n). Using the £"„ to reconstruct the tunnel- 
ing spectra, we find qualitative agreement with measured 
spectra [Q , and show that the i/-induced first-order tran- 
sition to the paramagnetic normal state observed for thin 
films can be softened in ultrasmall grains. 

Experimental Results. — Our goal is to understand in 
detail the i/-dependence of the measured discrete tun- 
neling spectrum (see Fig. 3 of [^) of an ultrasmall Al 
grain, coupled via tunnel barriers to one gate and two 
lead electrodes to form a nm-scale transistor. Each line 
in the spectrum corresponds to the i/-dependent energy 



Enj-E^r^^ + iE^-E^^^) needed for some rate-limiting 
electron tunneling process \ni)N±i |'^/)Ar off or onto 
the grain, where |n)7v is an eigenstate (with eigenenergy 
En + E^) of the iV-electron island with charging energy 
E^ . Since the change in charging energy 5Ec depends 
on the adjustable gate voltage Vg, so does the odd-even 
ground state energy difference {E^""^ - E^ ) (so that the 
BCS gap for the ground state of an odd grain can not be 
measured absolutely). Therefore only the energy differ- 
ences E^. — E^ I between same-A^ final states of tran- 
sitions with the same initial \ni)N±i-i he. the spacing 
between lines of a given spectrum, are physically signifi- 
cant. They give the grain's fixed- A^ eigenspectrum. By 
appropriately adjusting Vg, both even and odd spectra 
(A^ = 2m + p, with p — (0,1) for (e, o) parity) were mea- 
sured, and non-equilibrium effects minimized. 

The presence (absence) of a clear spectroscopic gap 
2il> d between the lowest two lines of the odd-to-even 
(even-to-odd) measured spectra (Fig. 3(a,b) of ||^) re- 
veals the presence of pairing correlations [|l|j4|: in even 
grains, all excited states involve at least two BCS quasi- 
particles and hence lie significantly above the ground 
state, whereas odd grains always have at least one quasi- 
particle and excitations need not overcome an extra gap. 

Pauli paramagnetism. — The measured levels' approx- 
imately linear behavior with H can be attributed to the 
electrons' Zeeman energies ±/i = ±^iJBgH; indeed, the 
differences between measured slopes of up- and down- 
moving lines correspond to (7-factors between 1.95 and 2. 
(Deviations from g = 2 probably result from spin-orbit 
scattering, known to be small but non-zero in thin Al 
films 0], but neglected below.) Thus, the i?-dependence 
of the spectroscopic gap i}{H) is almost entirely of Zee- 
man origin. Note that rt{H) must be distinguished from 
the BCS pairing parameter A(i?). In contrast to bulk 
samples, in ultrasmall grains the suppression of A(ff) 
through orbital diamagnetism is very weak, just as in 
thin films in parallel fields |^] : The flux through the grain 
(whose radius 5nm <C the penetration length of 50nm) 
is of order 5% of a flux quantum (pQ at H ~ 7T, i.e. too 
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small to significantly affect the orbital motion of electrons 
between reflections off the grain boundaries. Slight devi- 
ations from iJ-linearity observed in some larger grains ||] 
probably reflect the onset of orbital diamagnetism (giving 
corrections to eigenenergies of order ~ hvFr^{H/(j}o)^, [||). 
But for the spectra of interest here, they are much smaller 
than Zeeman effects and hence will be neglected. 

Now, Clogston and Chandrasekhar (CC) |^ argued 
that in the absence of orbital diamagnetism, supercon- 
ductivity will be destroyed by Pauli paramagnetism: Let 
|s) denote the ground state of the spin-s sector (s = 
J + p/2 with J an integer) of the (2m -I- p)-electron 
Hilbert space jl^, with exact eigenenergy Es{h,d) = 
Es{0,d) — 2sh ("spin" simply means J2^z with Sz = 
±^). CC pointed out that a ground state transition 
will occur from \p/2) to some normal state |s) when 
Ep/2(h,d) = Eg(h,d) at some sufficiently large field h. 
For d < A (as in thin films), Ep/2{0,d) = -A'^/2d 
and Es{0,d) — (s^ — p/4)d, with s ~ h/d [to ensure 
dsEs{h,d) = 0]. Thus, CC predicted a first-order tran- 
sition at a critical field hcc = A/-\/2, the new spin 
s — K/^/2d being macroscopically large. In tunneling 
measurements into thin (5nm) Al films (A = 0.38meV 
and Hqc — 4.7T) this first-order transition was observed 
as a jump in the tunneling threshold from A — /icc to at 
hcc- In contrast, the measured energy levels for ultra- 
small grains evolve continuously with h, showing kinks 
but no jumps. We suggest that this reflects a "soften- 
ing" of the transition that occurs when (i~ A, for which 
s should become of order one. We shall show this ex- 
plicitly by performing model calculations of the energies 
Eg{h, d) of BCS-like pair-correlated variational states \s). 

The model. — We adopt the reduced BCS Hamiltonian 
used in Ref. |Q] with an additional Zeeman term: 

H = J2i£j + ^hy^jaCja - \dY^c]^c]_Cj.-Cj,+ . (1) 

c|_|_ create free time- reversed states \j, ±), whose energies 
Sj = jd+Ef} — fi, measured relative to the chemical poten- 
tial /Lt, are taken uniformly spaced for simplicity (though 
this is not essential fl^). j = labels the flrst level 
whose occupation in the T — Fermi sea is not 2 but p, 
and \F) = 117=-™ c]_|_c]_ |Vac) is the even Fermi sea. 

Variational approach. — Since in the experiments T — 
50mK <C d, A, we set T — 0. For the spin-s ground state 
|s), we make a generalized BCS-like variational Ansatz 
| p2t , which pair-correlates all time-reversed states, except 
for 2s unpaired spin-up electrons placed as close as pos- 
sible to the Fermi surface, to minimize the kinetic energy 
cost of having more spin ups than downs: 

i^)= n svri("f^+^f^4-f4-)ivac). (2) 

]=~j j 

The prime over products (and over sums below) indicates 
exclusion of the singly occupied states j = — J to J+p— 1 



(for which u^^\v''^'> are not defined). Since (s|s') = Sss', 
the variational parameters v^'^'^ and uj'*'' must be found 
independently for each s (hence the superscript s), by 
minimizing the variational "eigenenergies" 

J+p-i 

8s{h,d) EE {s\H\s) = -2hs+ ^ Sj (3) 

j j j 

which we use to approximate the exact eigenenergies 
Es{n, d). The w^-term, not extensive and hence neglected 
in the bulk case, is non-negligible, but not dominant ei- 
ther. Solving dSs/dVj^^ — and u'^+v'^ = 1 simultaneously 
yields (wf '')^= {l-Sj /[e'^+Al]^/^)/2, with A^ determined 
by the generalized gap equation A^ — ^d^j and 

e^^Sj-Xdiv^f^)^. (To simplify our calculations, we used 
Ej =ef, this changes £s, which is stationary under small 
changes in v^^\ only to order A^. Then (TV) = 2m -I- p 
fixes the chemical potential to p = eQ — 5pfld/2.) 

The pairing parameter A^. — In the variational ap- 
proach. As is merely an auxiliary quantity in terms of 
which the v^j^^ and hence Eg are parameterized, and cer- 
tainly not directly measurable. It does serve as a mea- 
sure of the pairing correlations present in |s), though: if 
As — 0, |s) reduces to the paramagnetic spin-s ground 

state |s)o = n/=o~^ c]+ n7=-J '^J-l'^)' ^^'^ ^^'^ "corre- 
lation energy" \s\H\s) — o(s|77|s)o vanishes (see Fig. |l|b). 
The gap equation for As(d) is /i-independent (since the 
/i-dependence of Eg is so trivial), and differs from the 
standard bulk T = case (for which s = p/2, d<C A) 
through both the discreteness and the s-dependent re- 
striction on the sum, which respectively cause the d- 
and s-depcndence of As((i). Its numerical solution (see 
Fig. I^a) shows that As(d) decreases to zero as d is in- 
creased, because the kinetic energy cost of pairing corre- 
lations, which shift electron occupation probability from 
below to above ep, grows with increasing d More- 
over, As((i) decreases rapidly with increasing s at fixed d 
(reaching zero roughly at d = A/2s, as can be shown an- 
alytically). This s-dependence of Ag, a generalization of 
the even-odd effect (namely A1/2 < Ao) found in is so 
strong because for d ~ A there are only a few Sj 's within 
A of (where pairing correlations are strongest), so 
that increasing s and hence the number of unpaired elec- 
trons in this regime dramatically reduces the strength of 
pairing correlations. Evidently, for c? ~ A the standard 
mean-field description of superconductivity in terms of 
only a single pairing parameter A is no longer sufficient. 

Conceivably the £s , which are upper bounds on the ex- 
act spin-s ground state energies Eg, can be lowered by us- 
ing better variational wavefunctions that sample a larger 
portion of the spin-s Hilbert space, i.e. by "including 
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fluctuations" about the variational state |s). But in the 
present context this would hardly be worthwhile, since 
the Eg also depend quite sensitively on the unknown in- 
put energies {e^} [0, and here we merely seek a quali- 
tative understanding of the measured tunneling spectra. 
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of §). (b) The variational energies (s\H\s) - {F\H\F), 
plotted as functions of d at magnetic field h — (solid 
lines, ~ —A^/2d as d ^ 0) smoothly approach the ener- 
gies o(s|-H'|s>o -{F\H\F) = (s^ -p/4)d (dotted lines) of the 
uncorrelated states |s)o, since As{d) with increasing d. 
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FIG. 2. d-dependence of the critical fields 
induces, at fixed electron number A'', a ground state transition 
from \s) to |s'), labeled by (s,s') (in contrast we label the 
N-changing tunneling transitions in Fig. ^ by \n) W)). 
The lower curve gives the jump at ho^s{d) predicted for the 
lowest line of the e o tunneling spectra of Fig. ^. 

Critical Fields. — Having obtained As{d), we find the 
energies £s{h, d) numerically from Eq. (^). For two same- 
parity spins s' = s + J, £s' drops below £s at the field: 



,{d) = [£s'{0,d}~£s{0,d)]/{2s' ~2s) 



(4) 



Fig. |2| shows several ft-s,s' as functions of d. For given 
p, let their lower envelope be denoted by ^p/2,s('^) — 
min[ft,p/2,s'('^)]- This gives the "critical field" at which, if 
h is increased from at fixed d and p, the first change 
of ground state occurs from \p/2) to a new ground state 
\s{d)), whose spin depends on d. Numerically we find 
hp/^Ad ^ 0) = A/V2, which is CC's bulk resuh. We 
also find that for any d, \s{d)) always has As{d) — 0. 
Thus, the first ground state transition is always directly 
into a normal paramagnetic state with no pairing correla- 
tions. As anticipated above, the change in spin, s—p/2, 
macroscopically large at small d/A, decreases with in- 
creasing d. In this sense the first-order transition ob- 
served in thin films is "softened" for ultrasmall grains. 
s — p/2 reaches unity when the correlation energy be- 



comes smaller than the Zeeman energy gained by flip- 
ping a single spin: for p = (or 1) we find s = 1 (or 3/2) 
when d/A > OAT (or 0.32). This regime displays "min- 
imal superconductivity" : correlations are strong enough 
to cause a measurable gap, yet so weak that breaking a 
single pair destroys them. 

Tunneling Spectra. — In analogy to |s), one can also 
define excited states \s,n) in the spin-s sector of Hilbert 
space by placing the unpaired spins further away from 
the Fermi level, and variationally calculate their energies 
£s,nih,d) (writing £s = Ss.a)- It is easy now to recon- 
struct the expected tunneling spectra as function of h at 
fixed d, by finding the energy cost for all |si,ni)jv±i 
\sf,nf)N single-electron tunneling transitions. Since 
these satisfy the selection rule Sf — Si = ±1/2, only slopes 
of ±1 can occur. Neglecting non-equilibrium effects [§,D, 
we always take the ground state of a given spin-s sector 
as the initial state and denote it by \si{h,d),Q). The 
appropriate Si{h, d) must be determined from Fig. 

Whenever h passes through one of the critical fields 
hg. g', the current ground state of the iV ± 1 electron 
Hilbert space changes to |sj). This produces a discon- 
tinuity in the lowest line of the tunneling spectrum if 
Sf~s[ now violates the selection rule, or else a kink 
if only its sign changes relative to Sf — Si. According 
to Fig. H, depending on d one can distinguish different 
regimes I, II, III, . . . , in each of which the various Si 
to s[ ground state changes occur in a different order, 
leading to different magnitudes and positions of jumps 
in the tunneling spectra. In regime I, where the order 
of occurrence of ground state changes with increasing h 
is (0, 1), (i, |), (1, 2), (|, |), . . ., there are no discontinu- 
ities in the evolution of the lowest line [see Fig. ^(a,b)]. 
E.g. for the e ^ o spectrum, the lowest |0) |l/2) 
line changes continuously to |1) |l/2) at /io,i, since 
Sf — s[ — ^1/2. However, in all other regimes, s > 1 for 
the first ground state transition at /io,s, implying a jump 
in all e ^ o lines. The jump's magnitude for the lowest 
e — > o line, shown as function of d by the lower line in 
Fig. |, starts at d = from the CC value A(l - 1/V2) 
measured for thin Al films [Q, and decreases to (non- 
monotonically, due to the discrete spectrum), again illus- 
trating the softening of the transition. 

The absence of observable jumps in the measured low- 
est lines can be explained by assuming the grain to lie 
in the "minimal superconductivity" regime I |l3| . In- 
deed, the overall evolution of the lowest lines of Fig. 3 
of [^ qualitatively agrees quite well with those predicted 
for regime I [Fig. ||(a,b)], in which the correlation en- 
ergy is non-zero only for s — 0, 1/2. The prediction of 
ho^i ~ 0.95/icc is compatible with the experimentally 
observed value of Hq i — 4T, which is about 85% of 
Hqc = 4.7T for thin films [Q|. The jumps in higher lines 
(e.g. in Fig. |^(b) at are due to correlations left in 

exited states \s,n) (As^„ > A^^o since the unpaired elec- 
trons are further away from ep)- Experimentally, these 
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jumps have not been observed. This may be because 
up-moving resonances lose amphtude and are difficult to 
follow with increasing /i ^ , or because the widths of the 
excited resonances (« 0.13A) limit energy resolution 
More than qualitative agreement between theory and ex- 
periment can not be expected, for we assumed constant 
level spacing, neglected non-equilibrium effects, and the 
tunneling matrix elements are unknown. 




Magnetic Field /i/A 
FIG. 3. Tunneling spectra predicted for an ultrasmall su- 
perconducting grain for d/A = 0.67 and 0.34. Some lines are 
labelled by the tunneling-induced change s — + s' in the grain's 
spin. For clarity, not all higher lines are shown. Dashed lines 
indicate some of the critical fields hg^gi. 

Our theory predicts that for somewhat larger grains 
(in regimes II or higher) the tunneling spectra should 
show jumps even in the lowest line. It remains to be in- 
vestigated, though, whether orbital diamagnetic effects, 
which rapidly increase with increasing grain size (~ r^), 
would not smooth out such jumps. 

Non-time-reversed pairing. — By using the reduced 
BCS Hamihonian in Eq. (|l|), we neglected interaction 
terms -dj^ij^'j' Hh3,i' between non- 
time-reversed pairs cj^cj_, following Anderson's argu- 
ment B that for a short-ranged interaction, the matrix 



elements involving time-reversed states 



are much 



larger than all others, since their orbital wave-functions 
interfere constructively |14|. Interestingly, the experi- 
mental results provide strikingly direct support for the 
correctness of purely time- reversed pairing at /i = 0: if 
the X{j -\- 6j, j, j' -\- 6j, j') were all roughly equal to A for a 
finite range of (instead of being negligible for 5j 
as assumed in -ffred), then for 2s < 5j one could construct 
spin-s states with manifestly lower energy than |s), by 



pair-correlating non-time-reversed states, with the 2s un- 
correlated electrons at the band's bottom, where having 
them uncorrelated costs hardly any correlation energy: 



-m+2s-l 



n n ( 



t t 



|Vac) 



However, since it then costs no correlation energy to in- 
crease s (as can be checked explicitly), there would be 
no large threshold for h to induce ground state changes, 
and the change \p/2) to \p/2-\- 1) would occur at roughly 
h d (as in a normal paramagnet), contradicting the 
experimental finding that the first kinks only occur after 



a sizeable threshold h. 



p/2,. 



A/V2. 



Conclusions. — The dominant mechanism for destroy- 
ing pairing correlations in ultrasmall grains is Pauli para- 
magnetism. Calculating the variational eigenenergies of 
the lowest-lying eigenstates by a generalized variational 
BCS method, we have shown that decreasing grain size 
softens the first-order transition observed for thin films 
by reducing the number of spins flipped from being ma- 
croscopically large to being of order one for d ~ A. Our 
approach qualitatively reproduces the measured tunnel- 
ing thresholds, explaining why they do not show jumps. 
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